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^ ■ Abstract 



In this paper we prove a local-in-time existence theorem for an initial-boundary 
value problem related to a model of temperature-dependent phase segregation that 
generalizes the standard Allen-Cahn's model. The problem is ruled by a system of 
two differential equations, one partial the other ordinary, interpreted as balances, 
respectively, of microforces and of microenergy, complemented by a transcenden- 
tal condition on the three unknowns, that are: the order parameter entering the 
standard A-C equation, the chemical potential, and the absolute temperature. The 
results obtained in our recent paper [3] dealing with the isothermal case serve as 
a starting point for our existence proof, which relies on a fixed-point argument 
involving the Tychonoff-Schauder theorem. 
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1 Introduction 

This paper is a sequel and a generalization of our article [3], where we studied a nonlin- 
ear evolution system of the AUen-Cahn (A-C) type, intended to provide a mathematical 
description of the phenomenology of phase segregation by atom rearrangement in crys- 
talline materials, in the absence of diffusion. Our present generalization consists in taking 
thermal effects into account. To help understanding the underlying physics, temperature- 
independent mathematical models of A-C type and their proposed generalization are 
briefly discussed in the next section, where we recapitulate Gurtin's derivation of the 
standard A-C equation as well as the derivation of the nonstandard A-C system analyzed 
in [3], and where we sketch the main traits of our present temperature-dependent model. 
In Section 3, we formulate carefully the corresponding mathematical problem and we state 
a local-in-time existence result, that we prove in our last section. 



2 Phase segregation models of A-C type 

2.1 The temperature-independent model of |[3| 

In [5], we consider a nonlinear evolution system consisting of the partial differential equa- 
tion: 

Kdtp-Ap + f'{p) = fi (2.1) 

and of the ordinary differential equation: 

dt{-fi^p)=fi{^{dtpf + a), (2.2) 

complemented with the homogeneous Neumann boundary condition: 

dnP = on the body's boundary (2.3) 

(here dn denotes the outward normal derivative) and with the initial conditions: 

p\t=o = Po bounded away from , /i|t=o = /^o > . (2.4) 

The parabolic PDE (12.11) and the first-order-in-time ODE (12. 2p are interpreted as 
balances, respectively, of microforces and of microenergy. They are to be solved for the 
order-parameter field p = p{x,t) G [0,1], interpreted as the scaled volumetric density 
of one of two coexisting phases, and for the chemical potential field //. Moreover, k > 
is a mobility coefficient, / is a double-well potential confined in (0, 1) and singular at 
endpoints, and a = a{x, t) denotes a given source term. The microentropy field t] = —fj^^p 
cannot exceed the level from below, so that the corresponding prescribed initial field 

?7|t=o = '7o, (2.5) 

with ?7o = — /igpoj is nonpositive- valued. Taking fj, = in (12.11) yields the standard 
AUen-Cahn equation, which is intended to describe evolutionary processes in a two-phase 
material body, phase segregation included. 



2.2 Gurtin's derivation of the A-C equation 

The derivation of the A-C equation proposed by Gurtin [6] (see also [5] and [7] for similar 
derivations and for discussions of related models) is based on a balance of contact and 
distance microforces: 

div,C + 7r + 7 = (2.6) 

and on a dissipation inequality restricting the free-energy growth: 

dt^<w, w:=-7rdtp + ^-V{dtp), (2.7) 

where the distance microforce is split in an internal part it and an external part 7, the 
contact microforce is specified by the microscopic stress vector ^, and w is the (distance 
plus contact) internal microworking. Requesting the Coleman-Noll (C-N) compatibility 
[1] of the constitutive choices: 

7r = 7r(p,Vp,9ip), ^=i(p,Vp,9ip), and ^ = V^(p, Vp) = /(p) + ^|Vpp , (2.8) 

with the dissipation inequality (12. 7p yields: 

7r(p, Vp, dtp) = -f'{p) - k{p, Vp, dtp)dtp and |(p, Vp, dtp) = V p. 

Under the further assumptions that fi;(p, Vp, dtp) = n, a. positive constant, and that 7 = 0, 
the microforce balance (12. 6p takes the form of the standard AUen-Cahn equation 

Kdtp-Ap + f'ip) = 0. (2.9) 



2.3 The derivation of the A-C system of [3 



In [3], on adopting an approach put forward by one of us in [8], we deal with a modified 
version of Gurtin's derivation, in which the dissipation inequality (12. 7p is dropped and 
the microforce balance (12. 6p is coupled both with the microenergy balance 

dte = e + w, e:=—dwh + a, (2.10) 

and the microentropy imbalance 

dtT] >— div h + a, h := ph, a := pa (2-11) 

(here a is the external source of energy per unit volume). With a view toward modeling 
phase-segregation, we postulate that the microentropy inflow {h,a) be proportional to the 
microenergy inflow {h, a) through the chemical potential p, a positive field; consistently, 
we define the free energy to be: 

^:=e-p-V (2.12) 

with the chemical potential playing the same role as the coldness ■&~^ in the deduction of 
the heat equation. Combination of (I2.10p - (l2.12p leads to the inequality: 

dt^ < -T] dt{p~^) +p~^h-Vp-ndtp + $- V{dtp), (2.13) 



which replaces (12 .yp as a restriction on constitutive choices. We assume that, in addition 
to the independent variables p, Vp and dtp, the constitutive mappings delivering vr, ^, 77, 
and h depend also on /i; moreover, we choose: 

^ = ^(p, Vp, p) = -/ip + /(p) + \\^P?. (2.14) 

To impose the C-N compatibility of these assumptions with f l2.13p makes sense, because 
we have at our disposal two independent controls 7 and a to guarantee the free linear 
continuation in time of any given process t H- (p, p)(t) at any fixed space point. We find: 

7r(p, Vp, dtp, p)= p - /'(p) - k(p, Vp, dtp)dtp, |(p, Vp, (9tp, p) = Vp, 

r]{p,Vp,dtp,p)=-p'^p, h{p,Vp,dtp,p) = 0] (2.15) 

the last of these findings implies that the microenergy balance - in general, a PDE - 
becomes an ODE, a crucial mathematical simplification that we exploit in the following 
just as we did in [3]. Finally, under the additional assumptions that the mobility is 
a positive constant and the external distance microforce is null, the microforce balance 
(12. 6p and the energy balance (I2.1UJ) become, respectively, (12.11) and (12. 2p . 

2.4 Accounting for thermal effects 

As is well-known (see e.g. [5]), the classic heat equation can be arrived at by coupling the 
energy balance 

dte = -divh (2.16) 

and the entropy imbalance 

dtr]>-diYh, h = d^^h, (2.17) 

with the following constitutive prescriptions: 

^ = e-^r/, ^ = ^(^) = -c^^(ln^-l), (2.18) 

with the absolute temperature field "d positive-valued and the specific heat c^ a positive 
number. To account for thermal effects on the phenomenology of phase segregation by 
atomic rearrangement, we compare the formats (l2.16p -( E.18p and (l2.10p -( E.14p . A way 
to match them, in the light of the relationships of temperature and chemical potential to 
entropy provided by statistical mechanics, is to assume that: (i) the microenergy balance 
keeps the form (I2.10p : (ii) the energy /entropy fluxes and the free energy have the mutually 
consistent forms 

h = {d~^p)h, i) = e-{'dp-^)7], 7/^ = V^(p,Vp,p,^). (2.19) 

The second assumption is the main element of novelty of this note. With that measure, 
the dissipation inequality that replaces for (I2.13P is: 

dti^ < -V dt{^p~^) + {^p~^)h ■ V(^-V) -^dtp + i- V{dtp), (2.20) 



where, in addition to fl2.19l) -:i for the free energy, the distance force, the microscopic stress, 
the entropy, and the microenergy influx, are assumed to depend on the hst of variables 
A = {p,Vp,/i,^;9tp,V^}0 

When C-N compatibility of the present constitutive prescriptions with fl2.20p is in- 
spected, a delicate modeling issue emerges: this time, we cannot count on as many con- 
trols as needed to guarantee the free local continuation in time of any given process 
t I—)- (p, p, ^)(t) at any fixed space point. We do get the counterparts of the first, second, 
and fourth of f l2.15p . namely, 

7? = -dptp - ndtp, $ = d^pi) = Vp, h = 0, (2.21) 

and we are left with the residual inequality: 

{d^i; - ^fx-^v)dtf^ + {d^i^ + f^-'v)dt^ < 0- (2.22) 

Now, if it were possible to choose both dtfi and dt'& arbitrarily, then (I2.22p would yield 
the double equality: 

ry = ^-^n^dpijj = -p d^ij. (2.23) 

This observation motivates our decision to complement the microforce and energy balances 
with another field equation, namely, the thermodynamic consistency condition: 

fxd^ip + ^d^ip = 0. (2.24) 

With this, (I2.22p can be written as 

i9{d^ip - i^iJ.-^r])dt{^-^iJ,) < 0, (2.25) 

and (I2.23P follows, provided the time rate of (?9~^p) can be chosen arbitrarily. 
Next, we specify the free energy density ([2719]) 3 as follows (cf. (I2.14p ): 

V' = ^(p,Vp,p,^) = -pp + (p(p,^) + ^|Vp|^ (2.26) 

with 

<^(p, ^) = /(p) - c,^(ln ^ - 1) - Co p(^ - ^J, Co > 0, (2.27) 

where the double-well potential and the purely caloric free energy are supplemented by 
a coupling term that is effective when and where the temperature differs from the char- 
acteristic temperature "dc- With this final constitutive choice, the consistency condition 
IK2^ reduces to: 

pp + co^p + c„^ In ^9 = 0; (2.28) 

moreover, the balance of microforces (12. 6|) becomes: 

Kdtp-Ap + f{p)-co'& = ii (2.29) 



^Notc that the last two variables in the list give way to incorporate in the model the dissipation 
mechanisms relative to, respectively, atom-rearrangement without diffusion and heat conduction. Indeed, 
it is clear that ()2.19|) covers both special cases when either temperature or chemical potential is a space- 
time constant. 



where, with shght abuse of notation, we have written f'{p) for {f'{p) + coi)c)- As to the 
microenergy balance fl2.10p . we find: 

dt{-^-'fi^p) = ^-V(^ + ^^{^tpy), (2.30) 

an equation to be compared with (12. 2p . Our current mathematical model regards pro- 
cesses of phase segregation by atomic re-arrangement in the presence of thermal effects 
as solutions of the system of equations f l2.28p . fl2.29p and f l2.30p . 



3 Mathematical formulation and results 

The A-C system we derived is more difficult to deal with than its temperature-independent 
version we tackled successfully in [3], let alone the standard A-C equation (12. 9p . In 
fact, in addition to a PDE and an ODE as in [3], we now have to take care also of the 
transcendental equation (I2.28p . Our strategy is to repeat, for as much as is possible, the 
procedure in [3]: accordingly, we discuss the ODE (I2.30p first together with the relative 
initial condition, then we pass to the PDE and the transcendental equation, together with 
the relative boundary and initial conditions. 

3.1 Preliminaries 

In order to carry out the first part of our program, we adopt a change of variable to give 
(I2.30p plus (12. 5p the form of a parametric initial- value problem. We begin by introducing 
the initial value "i^o of ^ and setting: 

We then have that 

/^ = VU~P, (3.1) 

whence the Cauchy problem: 

^5,(^-^0 + ^^^4^ V^ = 0' (^~'0\t=o = -Vo. (3.2) 

Next, we restrict attention to the class of processes such that 

and replace (13. 2 p by the simpler problem: 

vP 

parameterized on both the space variable x and the field p{x, ■). Although simpler, this 
Cauchy problem is by no means trivial, because it can exhibit the Peano phenomenon and 
have infinitely many solutions; just as we did in [3], we pick a suitably defined maximal 
solution ^ (or y/^), having the important property to stay positive as long as is possible. 
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Remark 3.1. It remains to be seen whether the class we restrict attention to does include 
interesting phase-segregation processes. At this time, we cannot do any better than 
planning to check if this is the case by running numerical simulations. 

To complete our program, we note that, with fl3.ip . f l2.29p and f l2.28p become, respec- 
tively, 

K dtp - Ap + f'{p) - co'd = ^/Ifp , (3.4) 

and 

A(p,^) :=Cop^ + c„^ln^ = -V^pC", (3.5) 

an integro-dijjerential system for p and ?9, where ^/^ is implicitly defined in terms of p 
as the maximal solution to (13. Sp . This system is to be supplemented with the boundary 
condition ( 12. 3p . the initial condition for p in ( 12. 4p . and a compatible initial condition for ^. 

3.2 Results 

In view of the above discussion, we look for suitably smooth triplets of time-dependent 
fields {p,^,'d) over a regular region Q with boundary F, such that: 

< p < 1, ^ > 0, and ^ > 0; (3.6) 

atp-Ap + f(p)-Co^ = VeTp, with dnP = onT; (3.7) 

a.e+^^Ve-O; (3.8) 

A(p,^) = -v^; (3.9) 

p(0)=po, e(0) = eo, and ^(0) = ^o ; (3.10) 

and that, moreover, 

^ is maximal among the ^'s satisfying (13. 8p and the second of (I3.10p . (3-11) 

The problem's structure is the same as in [3], apart for the modifications due to the 
presence of the temperature variable ■&. Two items deserve a supplement of discussion. 
The first is that, just as in [3], we assume that 

0</ = /i + /2, where /i, /a : (0, 1) ^ M are C^.f^^nctions, (3.12) 

/i is convex, /I is bounded, lim/'(r) = — oo, and lim/'(r) = -|-oo, (3.13) 

r\0 r/"l 

with the constant CQi)c thought of as incorporated in /2(p). The second item has to do 
with the admissible choices of initial data: not only they must agree with (13. 6p . and hence 
satisfy 

0<Po<l, eo>0, ^o>0, (3.14) 

but also with (13. 9p . that is to say, they have to satisfy 

A(po,^o) = -v/p^- (3.15) 



To see what restrictions this last condition imphes on the choice of "do, it is convenient to 
study the function the function Ar : s H- A(r, s) = Cors + c^s Ins on (0, +00) for a given 
r G (0,1). Clearly, (i) A,, is strictly convex and tends to as s tends to 0; moreover, 
(ii) the equation \r{s) = has in (0, +00) a unique solution, that we denote by s(r); 
finally, (iii) A^ has a unique minimum point, denoted by s(r); in summary, for each fixed 
re (0,1), 

dX 

0<s(r)<s(r), A(r,s(r)) = 0, and ^^(r, s(r)) = 0. (3.16) 

OS 

A simple computation shows that 

s{r) = e'^'"*'', sir) = e'"*'' , and 

A(r, s(r)) = — c„e~^~'^*'', where c* := Co/c„ . (3-17) 

Therefore, a necessary condition for the existence of •&() is that 

VA)^ < Cye-^-''*''° a.e. in fi, (3.18) 



i.e., that sup C < 0, where ( := y/Po^o "~ c^e~^~'^*P°. If such a condition is satisfied, 
and if we want to solve (a.e. in Q) the equation X{pq{x),s) = for s, then uniqueness 
holds if ({x) = 0, and s{pq{x)) is the unique solution. Otherwise, if the strict inequality 
holds, then there are two solutions, the one in the interval {0,s{pq{x))) the other in 
(s(po(a;)),s(po(a:))). 

For existence of a local-in-time solution (p, ,^,^), a modest reinforcement of condition 
(I3.18P and a proper choice of ?9o suffice, namely, 

sup(v/Mo-c^e"^"'^*''") < and ^0 > s(po) a.e. in (]. (3.19) 

Under these assumptions, we can state the following result. 
Theorem 3.2. Assume that ( 13.12p -( l3.13p and ( 13.19P hold. Moreover, assume that 

a G L°°(fi X (0, +CX))), (ay G L°°(0, 00; L^(fi)), and po, ^0, ^0 e ^°°(fi); (3.20) 

PoeH'iQ), dnPo\r = 0, ApoeL°°{Q), inf po > 0, suppo < 1; (3.21) 

eo>0, ^/^,eH\n), A(^o,Po) = -Vp^- (3-22) 
Then, there exist T > and a triplet (p, ^, ^) satisfying: 

p G H\0, T; L\n)) n C°([0, T]; i/'(fi)); (3.23) 

peLP{0,T;W'^-P{n)) for each p< +00, S^p G L°°(Q), Q := f]x (0,T); (3.24) 

eGL°°(Q)nW^^'i(0,T;Li(f])), ^ G L°°(g), 9*^ G L°°(Q); (3.25) 

infp>0, supp<l, inf^>0, (3.26) 

and solving problem fl3.7|) - fl3.1ip . 

In the next section, this existence result is proved by a fixed-point argument. The 
method we use seems to be of some interest, because it relies on the application of the 
Tychonoff-Schauder theorem in a weak topology. 
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4 Proof of Theorem 13.2 



By (13.191) . we can start by choosing eo > such that 

V^Mo < c„e-^-^*^° - 2£o a.e. in n. (4.1) 

Moreover, as the result is local, we fix a reference final time T* > (e.g., T* = 1) and 
assume T < T* in the sequel. Our method is the following. By looking at (I3.7p - (l3.8p and 
to (13. 9p . separately, we construct two maps 

3^i:^^(p,0 and J^'-ip^O^^ 

with proper domains. Namely, the domain of 3^i is a convex set % depending on T and on 
a further parameter M, and the domain of 3^2 is the range Jioi 3^i. Then, we prove that a 
suitable choice of T and M ensures that the range of 3^2 is contained in %. This allows us 
to look for a fixed point of 3^2° 3^i- To this aim, we want to use the Tychonoff-Schauder 
theorem. For that reason, % will be endowed with some weak topology. We start to 
construct 3^i. The whole argument relies on the technique of [3] and the whole paper has 
to be revisited. This is done in the next steps. In particular, as in |3], we have to consider 
both the Cauchy problem obtained by coupling equation (13. 8p to the second (I3.10p and 
the Cauchy problem 

Q^^^\dtvl±^^ = and ^(0) = ^o (4.2) 



(i.e., p is replaced by v in (13. 8p ). where v satisfies 

V e D{^) ■= {v e H\0, T; L\n)) : v > 0, l/v e L^{Q)} (4.3) 

and possibly further conditions later on. Moreover, the map $ : -D($) — ?■ L°°(0,T; L^{Q)) 
is defined by $(w) = ^/^, where ^ is the maximal solution to the Cauchy problem (14. 2p . 

Remark 4.1. In the sequel, our notation is going to refiect dependences, if any, on such 
parameters as, say, T and M; however, possible dependences on problem data such as ^2, 
/, a, and the initial data, will not be displayed. 

The crucial constants and the maximum principle. We set: 

i!}^ := inf s(r) and 13* := sup s(r). (4.4) 

0<r<l 0<r<l 

A simple calculation yields: 

< 4 = e-(^+'=*) and ^ = 1; 
thus, we require that i) obeys the following double limitation: 

^, < ^ < r; (4.5) 
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in particular, the condition ini'd > (see fl3.26p ) will automatically hold true. Now, we 
notice that (13 .yp differs from the analogous one of [3] just for the presence of — cot? on the 
left-hand side. Therefore, even though such a term is space and time dependent, for a 
fixed •&, it can be seen as a part of the smooth perturbation /2 of the nonlinear term. Just 
by thinking of that, the construction of the crucial constants p*, p*, and ^* can be done 
exactly as in the quoted paper, provided that the definition of M2 is modified as follows 

M2:=sup{|/^(r)-Cos|: pG(0,l), s e {A,^)}. (4.6) 

Then, the analogues of Lemmas 4.1-4.3 hold in the present case, provided that we as- 
sume (14. 5p . Indeed, it suffices to read /2(p) — co'd in place of /2(p) in the proofs. 

The convex set. We define the set % as follows: 

X={^e L^{Q) : ^, dt'd eL°°{Q), A<^< ^*, \dt^\ < M}. (4.7) 

Clearly, % depends on T and on the a real parameter M, even though such a dependence 
is not stressed in the notation. Moreover, it is non-empty, convex, bounded, and closed. 

In the next steps, -(9 G X is given and we want to solve (I3.7p - (l3.8p and the first two 
equations (I3.10p for (p, by applying the procedures of [3]. 

LP estimates. We make a general observation. If p G (l,+oo) and some function z 
solves 

dtz-Az = ge LP{Q*) in Q* := Q x (0, T*) and dnZ = on the boundary, 
z{0) = Zq with Zq G L°°(f2), Azq G L°°{il), and (9„Zo = on the boundary, 

then, the following estimate holds 

||f^t^||LP(Q*) + W^W LP {0,T*;W^-P{n)) < C*):)! I|5'I|lp(Q*) + |ko||L°°(n) + ||A2;o||L°°(f7) j, 

where Cp depends on Q, T*, and p, only. Therefore, the same constant Cp yields 

||f^i2;||iP(Q) + ||2||iP(o,T;W2.p{o)) < C'pf ||5'I|lp(q) + ||2o||l=°(q) + ||A2;o|U°°(f7) j (4.8) 

for the solution z to the problem 

dfZ — Az = g E L^iQ) in Q, dnZ = on the boundary and z(0) = zq, (4.9) 

provided that T < T*. Now, note that z := p solves (14. 9 p with g '■= — f (p) + Cqi!) + \/^/p 
and that we are assuming {} E %; moreover, (I3.20p - p.21l) hold. Thus, by applying (14. Sp . 
we have: 

\\dtp\\LP{Q) + \\p\\LP{0,T;W^>P{n)) < Rp, (4.10) 

where Rp depends on the same parameters as in [3] and on -&*, but not on M or T (fol- 
lowing the rule laid down in Remark |4.H our notation stresses just the dependence on p). 
Furthermore, all this is true for the solution p coming from \/^ = ^{v), where v G -D($) 
satisfies p* < v < p*. 
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L°° estimate. As in [3], we differentiate with respect to time and see that u := dtp 
satisfies 

dtu — Am + u = F + o^dt'd and dnU = on the boundary, (4-11) 

m(0) = Apo - /'(po) + \/Up'o + co^o, (4.12) 

where F is as in |3], namely, 

F := dtp - f"{p)dtp -\^ p-^'^dtp -\x (|9,tf + -a) {vp)-''\ (4.13) 

where ip = $(t'), with v G -D($) such that p* < v < p* , and X is some characteristic 
function. Then, we fix the right value of g as in [3] , in order to get the desired estimate. 
We have 

IImIIl-cq) < C'^ll^ + CQdtd\\Li{Q) + ||m(0)||l-.(q) j 

where C does not depend on T (by the above general observation). As d E%, we have: 

\\F + Codtd\\L.iQ) < ||F|U,(Q) + Co(|fi|r)i/''M; 
by the use of the L^ estimates, we obtain: 

WdtPho^iQ) <Roo + C^MT^'\ (4.14) 

where R^o, C^c, and q are independent of M and T. 

The first map. At this point, for every ^ G X, we consider the following problem: in 

the set 



^={ip,0-- peH\0,T;L\n))nC'i[0,T];H\n)), ^eW''\0,T;L\n)) 

P*<P<P\ o<^<C 

\\dtp\\Lp{Q) + \\p\\LP{o,T;W^'P{n)) < Rp for every p E (1, +oo) 

Wdtph^iQ) <Roo + CooMT^/^}, (4.15) 

find a pair 

(p,0 satisfying ([32D, ^M, the first two of fl3lB . and (l3lTD (4.16) 

Indeed, the analogue of the map \1/ considered in [3] can be defined in the same way and 
actually maps its domain into itself also in the present case, because of our choice of the 
constants. Moreover, in performing the contraction estimate, just differences have to be 
considered and 'd disappears. Therefore, by setting: 

for ^ e%, Ji(^) is the solution (p,0 G ^ to problem (lilHjl . (4.17) 

we obtain a well-defined map 3^i : % ^- "Jl. We set: 

Jl is the range of 3^1. (4.18) 
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Such a subset of Ji depends on T and M and our next step is the definition of a map 
3^2 : 3i — ^ 3C that works as follows: given (p, ^ ^) the value 3^2ip,0 i^ ^ function {} 
satisfying (13. 9 P and the third of fl3.10p . However, before doing that, some more work is 
needed, since we are not sure that such a map is well-defined. It is not, indeed, unless T 
and M satisfy suitable constraints, and the next steps are devoted to find them. 

Higher regularity estimates. Let (p, ^) G 01. Then, (p, ^) solves f l4.16p for some 
•& & %, whence u := dtp solves fl4.1ip - p:.12p where F is the same as fl4.13p . On noting 
that f l3.2ip -( l3.22p imply m(0) G H^{Q), we have by the above L^ estimates and the general 
theory that 



\\9tP\\L^Q) + II f^tP i I L2{0,T;Hi(f^)) = \\9tU\\L^(Q) + \\u\\L2(^0,T:H\n)) 

< c{\\F + codtn^HQ) + WuminHn)) < C'{T, M), 



(4.19) 



where the form of the dependence of C'{T, M) on T and M is not important in the sequel. 

New a priori estimate. Let (p, ^) G 01. Then, (p, C,) solves (I4.16P for some t? G 3C. 
Therefore, as in [3], we have 



dtv^ 



-X 



\dtp\ 



a 



2v^ 



where X is some characteristic function, whence immediately 

1 



l^tV^II 



L°°(Q) 



< 



(ll^tPl 



+ m 



L°°(r2x(0,+oo)) 



)^ 



Hence, if c stands for different constants independent of M and T, we deduce: 



dtVp^ 
<cfl4 



^^+yp..ve 



^ v^ 11/^ II 



\dtpr 



i°°{Q) 






\dtv^\ 



L°°iQ) 



and we conclude that 



dt\^ <Ci(l + M2r2/«), 
where Ci is independent of T and M. 
First restriction on parameters. Observe that 

dt{v^ - c,e-^-'''')\ < \dt^/^\+c,c4dtp\. 
By accounting for (14.201) and (I4.14p . we deduce that 

Co 



where C2 is independent of T and M; this implies that 

||(V^-c,e-i-=*^) - (v/M^-c.e-i"'*''")|Uoo(Q) <C2T(l + Af2T2/'?). 
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(4.20) 



By (14. ip . we conclude that 

^/^ < c„e~^-^*P - eo a.e. in Q, (4.21) 

whenever M and T satisfy 

C2T(l + Af2T2/'') <eo. (4.22) 

Therefore, if T and M satisfy (g^S]), every (p,0 e 3^ fulfils (gJl])- 

The second map. Assume (p, ?9) G "51 and (I4.22p . By (14.211) . we have, in particular, 

that 

\fp\ < Cve~^~'^*^ a.e. in Q. 

Therefore, for a. a. (x,t) G Q, the equation A(p(x,t),s) = —\/p{x,t)S,{x,t) has two so- 
lution si and S2 which belong to (0, s(p(a;, t))) and (s(p(x, t)), s(p(x, t))), respectively. 
We term the latter ■d2{x^t) for a while and obtain a bounded function, thus a function 
■(92 G L'^iQ). Then, we can define 3^2 '■ "^ ^ L'^iQ) by setting: 3^2 (p, is such a function ■^2- 
Therefore, for (p, ^) G Jl, 

i!} = 3^2ip,0 means s{p) < 19 < s{p) and X{p,d) = —ypS, a.e. in Q. (4.23) 

As both p and ^ are continuous with respect to time (for a. a. x G f2) and the function 
{} := 3^2 (p,0 is always different from s{p), time continuity holds for -& as well, and it is 
clear that ^(0) = i^q. 

Although 3^2 is well-defined, we still have to find the restriction on T and M that 
ensures that the range of 3^2 is contained in %. 

Estimate from below. Let (p, ^ ^? set i!) := 3^2ip,0^ ^^^ assume (I4.22p . Then 
(I4.2ip holds. For almost every (x,t) G Q, we write the second-order Taylor expansion 
of the function s H- A(p(x,t),s) with center at s = s{p{x,t)) (hereafter, to lighten our 
notation, we refrain from writing {x,t)). We find s' G (s(p),?9) such that the following 
holds: 

r)\ 1 (9^ A 

- a/K = A(p, ^) = A(p, s{p)) + -(p, .(p)) (^ - s{p)) + 2 ^(P' ^') (^ - dp)? 



CtjC 



1-c.p , i^/,Q_ ,(^\\2 



lji^-s{p)f 



2s' 
in view of (13A6|1 and (13T7D . By (USD, we deduce that 



i^M- crn^^2 _ ^ ^-1-c.p 

whence 



^ i^ - s{p))' = c.e-'-''*' - Vp^ > Bo, 
2 s' 



[v — s[p)) > — s > — s[p) = — e *^ > — e ''^ . 

C^ Cy Cy Cy 

As i? > s(p), we conclude that 

^ - s{p) > 25o a.e. in Q (4.24) 

(with an obvious definition of 6o] note that 6o is independent of T and M). 
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Lemma 4.2. Assume that r G {p^, p*), s G {'d^,'d*), and 6 > 0. Then, 

dX 

s > s(r) + S implies t— (r, s) > Cy ln(l + 5/i!}*) . (4.25) 

OS 

Proof. Indeed, if s > s(r) + S, we have that 

(r, s) = Cor + Cy + Cyins = c^r + Cy + Cy In s(r) + Cyhi 



OS s\r) 

^"^f f W 1 ■^ 1 -^ ^ 1 sir) + 5 . , r, / x\ 

= Tr(^'^(^)) + Cy\xi—--- = Cyln—-- > Cyln --— = Cyln[l + d s{r)) 

OS s{r) s[r) s{r) ^ ' 

and the desired conclusion follows, because sir) <•{}*. D 

Estimate of the time derivative. Assume that (p, G 3^, ?9 := J^i^p,^, and 
that fl4.22p holds. On the one hand, we have that 

— {p^^)dt'& = -dt\/~pi- —{p,'&)dtp = -dt\/~pi - coddtp a.e. in Q. 

On the other hand, ( 14.24^ holds. Therefore, if we apply the previous lemma for 5 := 2So 

and set Lq := (c^, ln(l + 26o/'d*)) , we obtain: 



m\ < Lo{\dtVp^\+Comp\) < Lo{\^t^M\+Cor\^tp\). (4.26) 

Now, we account for ( 14.20p and conclude that 

\dt^\ < LoCi{l + M^T^/'^) + Loco^ (i?oo + CooMT'/") a.e. in Q. (4.27) 

At this point, we are ready to choose the constants T and M. 
Choice of the constants. Clearly, (14.271) implies |9(^| < M, whenever 

LoCi(l + M^T^/'^) + Loco^*(i?oo + CooMT^/^) < M. (4.28) 

Therefore, we choose, e.g., M = LqCi + LqCqi^* Roo + 1 and T small enough for both (14.221) 
and (14.281) to hold. With such a choice, the range of S^i is contained in Jl and 3^2 is a 
well-defined map whose range is contained in %. Hence, 3^ := 3^2 ° 3^i is a well-defined 
map from % into L'^{Q) that satisfies 3'(X) C %. 

Conclusion of the proof. As anticipated, we aim to an application of the Tychonoff- 
Schauder fixed point theorem. As far as the topology of % is concerned, we see IX as a 
subset of the topological vector space obtained by endowing L'^{Q) with the weak topology. 
Therefore, the convex set % is compact. So, the last point of our proof is the continuity 
of 5" with respect to the topology of %. To this aim, we observe that LF'{Q) with its 
strong topology is reflexive and separable and that % is bounded. Thus, the topology 
of % comes from a metric. In particular, 3" is continuous if and only if it is sequentially 
continuous. So, we pick 'dn,^ such that 

-^i G 3C for every n and ^n -^ ^ weakly in L'^i^Q) (4.29) 
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and pass to prove that the corresponding functions ^n '■= 3^{'&n) and ■& := S^ij)) satisfy 

4 -> ^ weakly in L\Q). (4.30) 

It suffices for us to show that every subsequence of {■/?„,} has a subsequence that converges 
to i} weakly in L'^[Q). In the sequel, to lighten our notation, we denote by the same symbol 
{{Vn}, say) both a sequence and all its subsequences. So, we start by any subsequence 
{'dn} of the given sequence and look for a subsequence of it satisfying fl4.30p . We set for 
convenience {pn,^n) '■= 3^i{^n)- Therefore, we have {pn,^n) E ^ C ^ and all the estimates 
given in definition fl4.15p hold true, as well as fl4.19p . In particular, we have 

Pn^ p weakly in W^'P{0, T; LP(fi)) n Lp(0, T; W'^'P{n)) for every p < +00 (4.31) 
dtpn -^ dtp weakly in H\0, T; L\n)) n L\0, T; H\n)) (4.32) 

for some p, at least for a subsequence. From (I4.31I) - (I4.32I) we also have that 

Pn — !■ p uniformly in Q and dtPn — ^ dtp strongly in L^{Q). (4.33) 

As p* < Pn < P* for every n, we infer the uniform convergence of f'{pn) to f'{p). Now, we 
consider the maximal solution ^ to the Cauchy problem obtained by coupling (13. 8p (with 
such a p) and the second (13.101) . We prove that (p, = '3^i{.'^)- By using [3l Lemma 4.8], 
we easily deduce that, a.e. in Vt and for every n and t G [0,T], 



WUt) 


-Vm 


Jo 


dtp, 


^is)' 


-\dtpis)\' ds + 


Jo 


\pn 


Jo 


dtpn{s) - 


dtp{s)\ds + i 


Jo 


\pn{s) - p{s)\ds, 







[s) - p{s)\ds 



where c stands for different constants independent of n. We infer that 

II V^ - V^I|l-(0,T;L1(O)) < c{\\dtPn - ^tp\\L^Q) + ||Pn " PIUhQ)' ) 



and, owing to (14.331) . we deduce that -y/^„ converges to -y/^ in L°°(0,T; L^(r2)), thus a.e. 
in Q at least for a subsequence. At this point, it is clear that {piC) satisfies problem 
fl3n)-f l3llD where we read ^ in place of d. Indeed, (ESD, the second f l310D . and f l3lTl) 
hold by definition of ^. On the other hand, it is straightforward to let n tend to infinity in 
the variational equation, or in an integrated version of it, and in the Cauchy conditions. 
We conclude that (p, = '3^i{'^)^ and the next step is the convergence of {On} to f? 
weakly in L'^{Q), at least for a subsequence. As a matter of fact, we can prove a strong- 
convergence result. Recall that the estimate from below (I4.24p holds for ^„. It follows that 

4 > s{Pn) + 2(5o = e''~'*P" + 26o = e-^-'*"" - e"^"^*'' + s(p) + 25o, 

whence 

"^11 > s(p) + 5q provided that |Pn — p| < ''705 

where r^o > is given by the uniform continuity of the exponential on every bounded 
interval (namely, r^r' G [p*^p*] and \r — r'\ < 770 imply that |e~^~'^*'' — e~^~'^*'' | < 6q). 
Moreover, the similar inequality 

s > s{r) + 60 
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holds true whenever \pn — p\ < Vo and 

min{p„, p} < r < max{p„,p} and min{'(9„,'(9} < s < maxj-i^^,-!?}. (4.34) 

On the other hand, the uniform convergence fl4.33p imphes that \pn — p\ < rjQ in Q ioi n 
large enough. Therefore, for such values of n, the following holds. We apply the Lagrange 
mean value theorem to A. For a.a. {x,t) G Q (once again, we omit writing {x,t) in the 
sequel) and suitable (space and time dependent) r„ and s„ satisfying (I4.34p . we have that 

VPn ^n - \/pi = KP^^) - KPn,^n) = -^(^n, S„) {p - Pn) + ■^(r„,S„) {^ - ^n) ■ 

Now, we observe that r„ G [p*, p*] and s„ > s{rn) + Sq since |p — r„| < rjQ. Hence, we can 
apply Lemma [4.21 and arrive at (14.251) with r = Vn, s = s„,, and 6 = 6q. This implies that 

ln(l + V^*) 14-^1 < Iv^P^- v^|+Co|pn-p|. 

As the right-hand side tends to a.e. in Q and is bounded in L°°{Q), we deduce that f^^ 
strongly converges to ?9 in L'^{Q). This concludes the proof. 
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